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A  COMPLEX  SYSTEMS  APPROACH  TO  METALLIC  FRACTURE 
SURFACE  CHARACTERIZATION 

Introduction 

The  physical  scale  of  observation  and  the  mathematical  formalism  of 
description  are  important  choices  for  the  modeling  of  physical  phenomena. 
For  continuous  media,  observations  at  different  scales  often  demonstrate 
distinctly  different  and  complex  features.  For  the  deformation,  damage  and 
fracture  of  solids,  observations  from  the  crystal  lattice  scale  to  the 
macroscopic  component  scale  define  the  domain  which  theoretical  and 
computational  modeling  address.  Within  these  boundaries,  quantitative 
relationships  between  any  intermediate  scales  are  valuable  aids  to  theoretical 
modeling. 

In  this  investigation,  the  geometry  of  the  large  scale  fracture  surface  is 
examined  in  terms  of  its  underlying  relationship  to  smaller  scale  geometric 
features,  quantified  at  approximately  the  continuum  scale.  The  measurement 
of  the  smaller  scale  fracture  surface  topography,  the  analytical  utilization  of 
these  measurements  and  the  role  of  computational  simulation  are  addressed  in 
the  course  of  demonstrating  large  scale  structure  on  the  fracture  surface.  This 
large  scale  structure  is  expressible  solely  on  the  basis  of  local  small  scale 
interactions. 

The  general  characteristics  of  the  large  scale  structure,  as  expressed  in  terms 
of  the  small  scale  orientations  of  the  surface  normals,  are  found  to  be 
relatively  insensitive  to  initial  conditions  applied  during  simulation.  This  is 
consistent  with  the  generally  reproducible  behavior  of  global  variables  such 
as  load,  displacement  and  crack  size  associated  with  multiple  material  specimen 
or  component  responses.  At  the  same  time,  the  uniqueness  of  small  scale 
feature  configurations  on  different  material  specimen  or  component  fracture 
surfaces  is  also  evident. 

This  discussion  is  organized  as  follows.  First,  factors  influencing  the  choice  of 
cellular  automata  and  video  imaging  as  analytical  and  experimental  methods 
in  this  investigation  are  outlined.  A  description  of  the  experimental 
arrangement  and  video  imaging  are  followed  by  a  section  on  the  calibration  of 
the  images.  A  discussion  of  the  cellular  automaton  formalism  is  followed  by  the 
image  processing  algorithm  which  derives  the  automaton  parameters  from  the 
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experimental  images  and  the  calibration  database.  Discussion  of  the  processed 
images  is  followed  by  the  generation  of  synthetic  images  and  comparison  to 
the  processed  images. 

Approach 

The  fracture  surfaces  of  metals  exhibit  diverse  features  at  different  scales  of 
observation.  Experimental  obervations  are  commonly  employed  across  the  full 
range  of  scales  from  the  crystal  lattice  to  some  portion  of  the  macroscopic 
component  size.  The  upper  limit  generally  depends  on  the  relative  proportion 
of  fracture  surface  to  component  size.  At  crystal  lattice  scales,  fracture  is 
viewed  in  terms  of  lattice  and  grain  boundary  interactions  with  irregularities 
in  these  structures,  e.g.  Garrison  and  Moody  [1],  At  microstructural  scales, 
fracture  is  described  in  terms  of  transgranular  or  intergranular  paths 

through  a  polycrystalline  field,  e.g.  Simmons,  Pao  and  Wei  [2]  and  Harvey, 
Sudarshan,  Louthan  and  Swanson  [3],  and  brittle  cleavage  or  ductile 
microvoids,  such  as  discussed  by  Thompson  [4]  and  Bassim,  Klassen,  Bayoumi 
and  Wilsdorf  [5],  Material  damage  evolution  at  or  near  the  the  continuum  scale, 
in  the  form  of  microvoids  and  microcracks,  as  observed  by  Chermant  and 
Coster  [6],  Bourcier,  Koss,  Smelser  and  Richmond  [7]  Becker  [8]  Becker, 

Smelser,  and  Richmond  [9]  and  Harvey  and  Jolles  [10]  for  different  stress  states 

and  materials,  is  followed  by  coalescence  at  the  mezzomechanical  scale 

discussed  by  Haritos,  Hager,  Amos,  Salkind  and  Wang  [11].  The  resulting 
macroscopic  crack  growth  exhibits  the  surface  features  generated  by  these 

cumulative  effects.  Recently,  fractal  geometry  has  been  examined  as  a  possible 

method  of  linking  the  irregular  geometries  of  fracture  across  spatial  scales  by 
Mandelbrot,  Passoja  and  Paullay  [12],  Underwood  and  Banerji  [13],  Pande, 
Richards,  Louat,  Dempsey  and  Schoeble  [14],  Pande,  Richards  and  Smith  [15] 
and  Dauskardt,  Haubensak  and  Ritchie  [16]. 

The  mathematics  used  to  classify  and  quantify  observations  at  each  scale 

should  reflect  both  the  geometric  features  and  analytical  forms  suggested  by 
the  observations.  The  periodicity  and  discrete  nature  of  atoms  in  a  body 
centered  cubic,  face  centered  cubic  or  hexagonal  close  packed  lattice  form  a 

concise  reference  on  which  to  identify  and  model  disturbances  associated  with 
substitutions,  interstitials,  vacancies,  dislocations  and  grain  boundaries. 
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Quantum  mechanical  electron  structure  calculations  are  explicitly  formulated 
to  reflect  these  underlying  reference  structures. 

The  presence  of  grain  boundaries  defines  another,  somewhat  larger  scale. 
Polygons  and  triangular  tesselations,  in  two  dimensions,  and  polyhedra  and 
tetrahedral  tesselations,  in  three  dimensions,  are  concise  topological 
structures  at  the  grain  scale.  Micromechanical  modeling  of  trans-  and 
intergranular  network  models  of  plastic  flow  and  fracture,  which  are 
explicitly  linked  to  these  intrinsic  geometries,  have  been  examined  by  Ostoya- 
Starzewski  [17,18]. 

At  continuum  scales,  associated  with  statistical  averages  of  the  local  anisotropy 
of  the  lattice  and  grain  scales,  continuous  and  differentiable  field  models  are 
traditionally  employed.  These  models  have  evolved  around  scalar  and  tensor 
field  descriptions  of  stress,  strain  and  energy  density  fields  intrinsic  to 
elasticity  and  phenomonological  theories  of  nonlinear  material  response. 
More  recently,  damage  fields  have  been  introduced  to  include  the 
phenomenological  effect  of  smaller  scale  features  in  a  continuum  formalism 
(e.g.  Kachanov  [19]).  Computational  models,  such  as  finite  element,  boundary 
element  and  finite  difference  techniques,  have  developed  and  incorporated 
these  concepts  to  facilitate  solution  of  two  and  three  dimensional  problems 
(e.g.  Curran,  Seaman  and  Shockey  [20]  and  Nemes  [21]).  Evolving  solutions  for 
nonlinear  boundary  conditions  such  as  fracture  surface  generation  and 
contact  area  size,  shape  and  friction  are  amenable  to  formulations  at  this  scale. 

Physical  measurements  of  metallic  fracture  surfaces  often  suggest  the 
evolution  of  smaller  scale  material  damage  into  fracture  surface  topographies 
composed  of  geometric  features  significantly  larger  than  the  continuum  scale. 
The  origins  of  these  features  lie  in  the  atomic  scale  lattice  inhomogeneities 

and  variations  in  grain  scale  sizes  and  shapes.  The  material  damage  evolution 

is  driven  by  the  stress,  strain  and  energy  density  fields  which  are  sensitive  to 

global  load  history.  Identical  locations  on  two  sets  of  fracture  surfaces, 
generated  from  two  specimens  of  identical  external  dimension  and  subjected  to 
the  same  load  history,  are  locally  unique.  Global  specimen  response,  however, 

is  often  quite  reproducible  from  specimen  to  specimen.  Therefore,  it  follows 
that  local  differences  in  material  damage  evolution  may  lead  to  differences  on 
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and  between  fracture  surfaces  which  may  not  translate  into  global 
differences. 

The  coupling  implied  by  such  local  and  global  responses  should  be  evident  on 
the  fracture  surface  at  a  scale  intermediate  between  the  continuum  and 
macroscopic  scales.  At  such  an  intermediate  scale  the  fracture  surface 
uniqueness  should  be  quantifiable  within  the  context  of  an  appropriate 

mathematical  description.  Experimental  measurements  to  support  such  a 
description  should  be  of  sufficient  spatial  resolution  to  delineate  local  relevant 
geometric  parameters  in  quantities  sufficient  to  facilitate  global  comparison 
in  a  statistical  sense.  Such  a  description  should  also  be  able  to  identify  surface 
equivalence. 

The  choice  of  a  compact  mathematical  description  to  parameterize  a  three 
dimensional  fracture  surface  is  influenced  by  the  nature  of  the  surface  and 
the  information  desired  from  the  surface.  Exact  characterization  of  each  local 

feature  on  a  surface  is  one  possible  objective  of  parameterization.  The 

identification  of  dominant  types  of  local  features  comprising  a  surface  is 
another.  A  method  of  equivalence  classification  between  surfaces  is  another 

possible  objective.  This  investigation  will  focus  on  the  latter  two  objectives. 

As  discussed  above,  the  physical  nature  of  the  fracture  process  suggests  that 

the  coupling  between  the  local  and  global  repsonses  takes  place  through  the 
local  interactions  between  the  material  deformation  and  damage  leading  to 
fracture.  An  appropriate  mathematical  form  to  incorporate  these  effects  is  a 

ztochastic  cellular  automaton.  Cellular  automata,  as  discussed  in  detail  below, 
account  for  local  interactions  and  operate  iteratively  to  generate  a  global 
configuration  as  a  set  of  discrete  site  values  on  a  mathematical  lattice.  Such  a 
representation  is  a  compact  and  efficient  method  of  description  and  facilitates 
computational  simulations  of  the  fracture  surfaces. 

The  underlying  restriction  of  cellular  automata  implementation  is  the  need  to 

limit  both  the  size  of  the  parameter  set  used  to  describe  the  lattice  sites  and  the 

size  of  the  site  neighborhood  used  to  account  for  local  interaction.  For  the 
purposes  of  this  investigation,  the  local  surface  inclination  was  expressed  as  a 
value  indicating  the  nature  of  orientation,  through  its  local  normal,  rather 
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than  its  specific  value.  Three  local  surface  inclination  classes  were  used:  null, 
principal  and  diagonal  (Figure  1).  All  were  defined  (i)  with  respect  to  a 
threshold  angle  parameter  measured  from  the  macroscopic  fracture  plane  (i.e. 
its  elevation)  and  (ii)  the  direction  of  crack  growth  (i.e.  its  azimuth). 

A  null  inclination  was  defined  for  any  iocal  surface  orientation  for  which  the 
angle  between  the  local  surface  normal  and  macroscopic  surface  normal  was 
less  than  the  threshold  angle.  A  principal  local  surface  inclination  was 

defined  for  orientations  exceeding  the  threshold  and  being  inclined  either 
along  or  perpendicular  to  the  direction  of  crack  growth.  A  total  of  four 
different  local  surface  inclinations  are  equivalent  and  fall  under  the 

principal  classification  by  this  definition.  A  diagonal  local  surface  inclination 
was  defined  as  exceeding  the  threshold  but  falling  between  the  principal 
inclinations.  Again,  a  total  of  four  inclinations  are  equivalent  and  classified  as 
diagonal  inclinations  by  this  definition. 

Given  these  approaches  and  definitions,  the  objective  of  the  this  investigation 
was  to  examine  the  global  and  local  relationship  between  local  mezzoscale 
features  of  a  fracture  surface.  Video  images  of  the  fracture  surface  were  used 
to  set  the  spatial  resolution  scale  and  acquire  intensity  data  from  which  the 
surface  topography  will  be  assessed.  The  local  surface  orientation  was 

characterized  from  the  intensities  of  the  video  images  by  transformation  to 
local  surface  inclination  parameter  set  at  the  physical  scale  corresponding  to 
pixel  resolution. 

The  processed  image  was  used  to  identify  stochastic  cellular  automata 
parameters.  The  use  of  stochastic  cellular  automata  was  motivated  by  the  issues 
discussed  above,  i.e.  local  fracture  surface  uniqueness  and  global  fracture 
surface  equivalence  consistent  with  global  specimen  response  reproducibility. 
It  also  facilitates  future  modeling  and  simulation  of  material  damage  and 

fracture  surface  evolution. 
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Video  imaging  configuration 


Evolving  computational  capabilities  simultaneously  influence  theoretical  and 
experimental  strategies  to  support  the  modeling  of  physical  processes.  The  use 

of  video  imaging  for  quantitative  fractography  is  an  area  where  the  ability  to 
acquire  and  store  visual  data  in  sufficient  quantity  exists  to  support  the  data 

requirements  of  computationally  intensive  modeling  paradigms.  To  address  the 
issues  posed  at  the  outset  of  this  investigation,  video  imaging  was  used  to 

acquire  and  store  black  and  white  images  of  an  HY-100  compact  fracture 
specimen  of  standard  geometry.  The  specimen  was  precracked  at  room 
temperature,  cooled  in  liquid  nitrogen  and  rapidly  fractured.  (It  should  be 
noted  that  a  fracture  surface  generated  by  any  means  was  amenable  to  the 
imaging  and  subsequent  analysis).  A  40  by  70  pixel  image,  covering  an  8.41  x 
10‘2  cm2  area  in  the  center  of  the  fracture  surface,  was  used  for  this 

investigation.  Each  pixel  corresponded  to  approximately  3.04  x  I  O'5  cm2  of 
projected  area  on  the  fracture  surface,  or  approximately  5.51  x  10‘3  cm  per 
pixel  side. 

Each  of  the  2800  pixel  intensities  composing  the  image  were  stored  as  8-bit 
values  providing  an  intensity  range  from  0  to  255.  These  images  were 
processed,  in  the  manner  discussed  below,  to  express  the  local  surface 

topography,  on  a  pixel  by  pixel  basis,  in  terms  of  local  surface  inclination 

classes.  The  local  surface  inclination  class,  as  discussed  above,  represented 
null,  principal  and  diagonal  inclinations. 

The  fracture  specimen  was  positioned  in  line  with  the  optical  axis  OP  of  the 
camera  (Figure  2).  The  plane  defined  by  the  fracture  surface  F,  when  the 
fracture  surface  normal  was  coincident  with  the  optical  axis  of  the  camera, 
was  defined  as  a  fixed  reference  plane  Rp  from  which  specimen  luiaiions  were 
set  for  the  calibration  images.  The  specimen  itself  could  be  rotated  about  either 
one  of  two  orthogonal  axes,  i.e.  the  N-S  and  the  E-W  axes.  Directions  will  be 
specified  by  angle  values  for  conciseness  and  common  directions,  i.e.  north 
(N),  south  (S),  east  (E)  and  west  (W),  for  case  of  discussion.  The  N  direction 

coincided  with  the  direction  of  crack  growth.  The  N-S  axis  definition  follows 

immediately.  The  E  and  W  directions  were  accordingly  transverse  to  the 
direction  of  crack  growth  and  define  the  E-W  axis.  Therefore,  the  axes  are 

6:0°.  90°  }  (I) 
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Schematic  diagram  of  specimen.  C'C'l)  camera  and  eight  illumination  directions 


e  =  {  e  :  E-W  ,  N-S  } 


(2) 


o  r 


and  lie  in  the  reference  plane. 

Two  sets  of  video  images  were  obtained  by  a  CCD  camera.  The  field  of  view  on 
the  fracture  surface  is  shown  in  Figure  3.  The  CCD  camera  optical  axis  and 
fracture  surface  reference  plane  normal  were  parallel.  The  first  set  of  images, 
set  C  =  {ce^y},  were  for  calibration  of  grey  scale  intensities  required  to  obtain 
local  surface  orientations.  The  second  set  of  images,  set  Q  =  {q<(>},  were  for 
processing  of  pixel  intensities  into  local  surface  orientations  and  assignment 
into  local  surface  inclination  classes.  The  subscript  $  refers  to  the  illumination 
source  and  belongs  to  the  set 


<t>  =  {  0  :  0°  ,  45°  ,  90° .  225°  .  270°  .315°  }  (3) 


o  r 


<t>  =  {  4> :  E  ,  NE  ,N . SW  ,  W  ,  NW  }  ( 4 ) 

The  subscript  y  refers  to  the  calibration  orientation  and  belongs  to  the  set 

V  =  {  y  :  -45.00°  ,  -33.75°  ,  -22.50° .  22.50°  .  33  75°  ,  45.00°  } 

(5) 

Illumination  of  the  fracture  surface  for  the  calibration  and  processing  images 
was  provided  by  one  of  eight  different  light  sources  in  the  illumination  set  L  = 
{ 1  (jj ) .  The  selection  of  the  particular  light  source  was  consistent  with  the 
requirements  of  the  calibration  and  processing  procedures  described  in  detail 
below.  The  light  sources  were  mounted  in  a  coplanar,  evenly  spaced  circular 
pattern  (Figure  2).  The  positions  of  the  light  sources  in  the  light  source  plane 
Rl  coincided  with  the  E,  NE,  N,  NW,  W,  SW,  S,  and  SE  directions  define  similarly 
to  those  on  the  fracture  surface  reference  plane  Rf-  The  normal  to  the  light 


9 


ig  '  Video  image  of  fracture  surface  Crack  growth  direction  is  from  bottom. 

toward  top  ol  image 
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source  plane  was  positioned  to  lie  on  the  optical  axis  O  of  the  camera.  The  light 
source  plane  was  placed  at  a  distance  away  from  the  fracture  surface  plane  so 
that  each  source  was  positioned  at  an  elevation  of  45°  with  reference  to  the 
point  defined  by  the  optical  axis  O  and  the  fracture  surface  reference  plane 
Rp.  Illumination  azimuth  was  therefore  controlled  by  source  selection  while 
the  elevation  was  fixed. 

Each  pixel  of  a  fracture  surface  image,  i.e.  cilte$y  or  qil,0,  was  the  image  of  a 
small,  local  section  of  the  fracture  surface.  The  i  and  t  variables  denote  the 
local  directions  parallel  and  perpendicular  to  the  crack  front.  These  coincide 
with  the  N-S  and  E-W  axes,  respectively.  The  calibration  pixel  intensities  were 
influenced  by  specimen  rotation  y  (Figure  4)  around  the  6  axes  and 
illumination  source  +  .  Preliminary  studies  of  the  fracture  surface 
demonstrated  its  local  pixel  scale  and  global  image  scale  diffusivity  toward 
incident  light.  This  was  used  for  calibration  and  processing  in  such  a  maimer 
that  no  assumptions  on  the  specific  character  of  the  diffusivity  were  required. 
Surfaces  of  a  more  specular  nature,  while  not  studied  in  this  investigation, 
could  also  be  amenable  to  study  in  their  natural  state  or  by  appropriate 
deposition  of  a  thin  diffusive  coating  which  would  not  alter  the  surface  at  the 
scales  of  interest  defined  by  pixel  resolution. 


Calibration  Procedure 

Calibration  of  the  fracture  surface,  through  the  image  set  C.  was  required  to 
establish  reference  intensities  with  which  to  compare  the  intensity  of  images 
for  processing  on  a  pixel  by  pixel  basis.  A  number  of  different  approaches 
have  been  used  to  address  the  characterization  of  surface  geometry  from  two- 
dimensional  images  of  reflected  surface  intensity,  as  discussed  in  the 

collection  of  approaches  to  this  class  of  problem  edited  by  Horn  and  Brooks 
[22].  In  many  applications,  such  as  terrain  mapping  and  image  recognition, 
only  a  single  image  or  a  pair  of  stereo  images  are  available.  The  object  being 

viewed  is  generally  not  under  the  control  of  the  observer.  The  approach  in 

this  investigation  drew  on  some  of  the  concepts  employed  in  the  study  of  these 
problems  while  taking  advantage  of  the  ability  to  acquire  images  of  the 
fracture  surface  in  different  orientations.  Fracture  surface  orientation  was 
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specimen  rotation 


varied  by  rotation  of  the  specimen  with  respect  to  the  N-S  axis  and  by  rotation 
with  respect  to  the  E-W  axis.  Different  conceptual  approaches  to  the 
calibration  problem  are  reviewed  along  with  the  actual  approach  implemented 
in  this  investigation. 

In  principle,  a  local  fracture  surface  calibration  would  compare  pixel 
intensities  against  reference  values  associated  with  the  pixel  itself.  These 
could  be  constructed  for  each  pixel  by  illuminating  the  fracture  surface  from 

different  directions  while  maintaining  a  constant  distance  and  orientation 
from  light  source  to  surface  to  ensure  constant  incident  intensity.  An 
intensity  calibration  over  both  illumination  azimuth  and  elevation  would  be 
obtained  for  each  pixel.  This  requires  the  camera  and  fracture  specimen  to 

move  as  a  system  if  the  light  sources  are  fixed,  or  the  light  sources  to  move  if 
the  camera  and  fracture  surface  system  are  fixed,  in  order  to  maintain 
correspondence  between  each  pixel  and  an  associated  local  region  of  the 
fracture  surface  through  multiple  image  acquisitions.  Practical  difficulties 
were  associated  with  being  able  to  move  either  the  video  camera  and  fracture 

specimen  as  a  unit  or  being  able  to  move  the  light  sources.  Although  this 
approach  was  not  used  in  this  investigation,  a  system  to  do  this  is  currently 
under  development. 

An  alternative  to  the  pixel  intensity  calibration  procedure  would  be  to  achieve 
a  global  surface  intensity  calibration.  Local  intensities  would  be  compared 
against  a  reference  values  reflecting  an  average  intensity  incorporating  the 
contributions  of  the  entire  region  under  consideration.  Conceptually,  such  a 

calibration  could  be  obtained  by  moving  the  camera  to  a  sufficiently  large 

distance  away  from  the  specimen  such  that  the  region  of  interest  on  the 

fracture  surface  would  be  imaged  in  only  one  pixel.  The  illumination  source 
could  again  be  positioned  relative  to  the  specimen  and  camera  to  cover  the 

required  range  of  azimuth  and  elevation.  The  difficulty  of  moving  the  camera 
to  a  sufficiently  far  distance  and  imaging  the  fracture  surface  region  of 
interest  on  a  single  pixel  are  obvious. 

A  practical  solution  to  the  calibration  problem  was  implemented  by  acquiring 
images  with  the  specimen  at  various  rotations  about  each  rotation  axes.  In  bis 

manner  both  the  camera  and  the  light  sources  were  stationary.  Processing  of 
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pixel  intensity  values  c; l,8<j>y  was  performed  for  each  pixel  in  each  image.  Each 
image  was  defocused  by  replacing  the  center  pixel  value  of  all  three  pixel  by 

three  pixel  subimages  by  the  average  of  the  nine  pixel  values.  This  is, 
essentially,  a  low  pass  spatial  filter.  Thus,  an  essentially  uniform  intensity 

distribution  over  the  fracture  surface  region  of  interest  was  achieved.  This 

process  was  analogous  to  the  far  field  calibration  described  above.  Direct 
global  averaging  of  the  pixels  to  a  value  C0$y  were  found  to  approximate  the 
successive  local  averaging  procedure.  On  this  basis,  direct  global  averaging 
was  adopted  for  purposes  of  surface  calibration. 

Calibration  images  were  obtained  for  specimen  rotations  about  the  E-W  and  N-S 
axes.  For  each  of  the  eight  illumination  directions  0  and  two  specimen  rotation 
axes  0,  nine  specimen  orientations  y  ranging  from  -45°  to  +45°  in  11.25° 
increments  were  illuminated  and  stored  for  a  total  of  144  calibration  images. 

The  projected  area  of  the  fracture  surface  region  of  interest  varied  since  the 
specimen  was  positioned  in  front  of  the  Fixed  camera.  The  image  processing 
region  was  modified  for  each  image  to  contain  only  the  region  of  interest  on 

the  fracture  surface. 

The  calibration  curves  obtained  from  this  global  averaging  procedure  are 

shown  in  Figure  5.  The  curves  are  plotted  for  a  constant  illumination 
direction,  with  respect  to  the  fracture  surface  plane  and  direction  of  crack 
growth,  from  measurements  taken  at  the  different  fracture  surface 
orientations  about  the  vertical  and  horizontal  axes.  Linear  interpolation  is 
used  to  link  together  the  nine  specimen  rotation  values  for  each  illumination 
direction. 

The  anticipated  monotonic  decline  of  the  illumination  intensity  ce<j>y  is 

apparent  as  the  fracture  surface  normal  vector  is  rotated  away  from  light 
source  direction.  An  example  is  the  E  illumination  calibration  curve  for  the 
the  specimen  rotation  about  the  vertical  N-S  axis.  Essentially  constant 

illumination  intensity  is  apparent  for  light  source  directions  and  specimen 
rotation  axes  defining  a  plane  perpendicular  to  the  plane  defined  by  the 
specimen  normal  vectors  as  the  plane  is  rotated.  The  N  illumination 
calibration  curve  for  the  vertical  specimen  rotation  axis  is  an  example.  The 
intermediate  illumination  directions,  such  as  NE,  exhibit  an  attenuated 
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Intensity  versus  specimen  orientation  calibration  curves  for  eight  light  sources. 
Specimen  rotation  about  (a)  vertical  axis. 


224  h  SPECIMEN  ROTATION  AXIS  -  HORIZONTAL 


16 


(Continued)  Intensity  versus  specimen  orientation  calibration  curves  for  eight 
light  sources.  Specimen  rotation  about  (b)  horizontal  axis. 


monotonic  decline  as  the  fracture  surface  is  rotated  away  from  the  light 
source  direction.  Symmetry  and  antisymmetry  of  magnitude  and  slope  is 
clearly  produced  between  appropriate  pairs  of  calibration  curves  as  expected 
for  the  proper  combinations  of  specimen  rotation  axis  and  illumination 
direction. 


Cellular  automata  theory 

Cellular  automata  were  discussed  by  von  Neumann  [23,24]  in  order  to  study  the 
complexity  and  means  by  which  objects,  assembled  from  simple  components, 

could  self-replicate  within  a  large  system.  The  discrete  nature  of  cellular 
automata  in  both  space  and  time  provides  a  large  number  of  identical  sites  in 
space  whose  time  evolution  is  governed  by  identical  local  neighborhoods  of 
adjacent  sites.  These  sites  are  usually  arranged  in  a  periodic  lattice.  Simple 
rules  govern  the  synchronous  evolution  of  each  possible  neighborhood  at  the 
current  time  to  a  new  site  value  at.  the  next  value  of  time.  The  value  of  each  site 
on  the  resulting  lattice  in  space  and  time  is  drawn  from  a  finite  set  of  site 

values.  The  outcome  for  each  possible  neighborhood  configuration  may  be 

unique,  in  the  case  of  deterministic  cellular  automata,  or  drawn  from  all 
possible  values  in  a  nonuniform  manner,  as  in  the  case  of  stochastic  cellular 
automata.  Large  scale  digital  computation  and  computer  graphics  has  made  the 
study  of  large,  synchronous  systems  more  accessible,  facilitating  interest  in 
modeling  paradigms  such  as  cellular  automata.  A  brief  review  of  cellular 

automata  follows  in  order  to  introduce  stochastic  cellular  automata  for  this 
investigation. 

Deterministic  cellular  automata  (e.g.  Wolfram  [25])  and  stochastic  cellular 
automata  (e.g.  Grinstein,  Jayaprakash  and  Y.  He  [26]  ,  Kinzel  [27],  Kaneko  and 
Akutsu  [28])  may  be  formulated  on  lattices  whose  site  values  form  a  finite  set. 
Define  xj{  to  be  the  site  value  at  position  i  and  time  t.  The  site  values  are 
defined  by  the  set  X  of  k  integers  as 

X  =  {  x  :  0  ,  1 . k-1  }  (6) 


so  that 
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(7) 


Xj*e  X 

Site  values  at  time  t+1  are  governed  by  an  automaton  rule  which  maps  the 
preceding  neighborhood  configuration  at  time  t  into  the  new  site  value  at  time 
t+1  (For  a  spatial  radius  r  equal  to  one,  the  neighborhood  is  shown  in  Figure  6). 
The  "radius"  of  the  spatial  neighborhood  is  defined  as  r. 

In  the  deterministic  case 

*it+1  =  F[  Xi.rl  ,  xi-r+il .  xjl  .....  Xi+r-i1  .xj+r1  ]  (8) 


o  r 


*it+1  =  F[  yit  ] 


(9) 


where  y,1  is  the  neighborhood  configuration  number,  i.e. 


Zr 


y,  =  2> 

fO 


2r-j  i 


(10) 


A  total  of  2r+l  sites  compose  each  site  neighborhood  at  time  t.  The 
configuration  of  the  x,1  neighborhood  governs  xit+1  through  the  application 
of  the  rule  F.  The  value  of  y,1  is  contained  in  a  set  Y  defined  by  the  size  k  of  the 
set  X  and  the  neighborhood  radius  r,  i.e. 

Y  =  (  y  :  0  ,  1  . k2r+1  }  (11) 

so  that 

y,le  Y  (12) 

The  neighborhood  configuration  number  is  a  useful  representation  when 
implementing  the  cel!",ar  automaton  rule  in  tabular  form.  The  table  required 
for  a  deterministic  cellular  automaton  would  be  a  list  k2r+1  elements  long,  i.e. 
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one  outcome  for  each  element  y  of  neighborhood  configuration  set  Y.  Tabular 
implementation  is  also  computationally  efficient  and  admits  all  possible  rules. 

The  exponential  growth  of  neighborhood  configurations  with  neighborhood 
radius  generally  constrains  the  size  of  the  set  X  for  statistical  and 
computational  reasons.  The  coarseness  of  the  set  X  and  the  neighborhood 
configuration  size  are  offset  by  the  large  number  of  rule  iterations  which 
govern  the  lattice  site  values.  A  site  value  depends  on  the  site  values  contained 
in  a  space-time  cone  which  converges  on  the  site  value  through  increasing 
time.  Simultaneously,  a  site  value  affects  any  site  in  a  space-time  cone  which 
expands  through  increasing  times. 

In  the  case  of  stochastic  cellular  automata,  the  site  value  set  X  and  finite  spatial 
neighborhood  are  identical  to  those  defined  for  deterministic  cellular 
automata.  The  configuration  y;1  formed  by  the  x,1  neighborhood,  however,  is 
associated  with  a  set  of  probabilities  py>x  governing  the  generation  of  the  site 
value  xil+1.  The  value  of  x,l+1  generated  from  the  neighborhood  configuration 
will  not  be  unique,  which  differs  from  deterministic  cellular  automata.  The 
stochastic  cellular  automata  rule  P  is 

Xit+1  =  P[  Xj.rl  ,  Xj.r+11 .  x;‘  .  Kj+r-11  .*i+rl  = 

Py.O  •  Py,l*  •••  •  PyJc-2  •  Py,k-1  ]  (13) 


o  r 


xjt+1  =  F[  yjl  ;  pyiX]  (14) 

The  probabilities  for  each  neighborhood  y  must  satisfy  the  requirement  that 

Py.O  +  Py.i  +  ■••+  Py,k-2  +  Py,k-I  =  ’•  (15) 

The  table  required  for  a  stochastic  cellular  automaton  would  be  an  array  of 
k  2r+ 1  rows  and  k  columns,  representing  the  probability  of  each  possible 
outcome  for  each  neighborhood  configuration  in  set  Y. 
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Stochastic  cellular  automata  for  fracture  surface  image  processing 


The  fracture  surface  is  physically  generated  in  a  sequential  manner  as  the 
crack  propagates.  The  propagation  direction  defines  a  direction  of  increasing 
time  and  a  complimentary  direction  of  constant  spatial  position.  The  two 
dimensional  spatial  domain  in  the  video  images  may,  therefore,  also  be 
interpreted  as  the  time  evolution  of  an  essentially  one-dimensional  linear 
portion  of  the  crack  front  as  it  propagated  into  the  solid  and  generated  the 
crack  surface.  Lines  of  constant  time  may  be  considered  parallel  to  the  crack 
front,  while  lines  of  constant  position  would  be  orthogonal  to  the  lines  of 
constant  time.  For  the  video  images  in  this  investigation,  constant  time  is  in 
the  horizontal  direction,  time  is  increasing  in  the  upward  direction  and 
constant  position  is  in  the  vertical  direction. 

Within  the  context  of  the  stochastic  cellular  automata  formalism  described 
above,  the  specific  form  used  for  the  fracture  surface  analysis  is  defined.  The 
local  surface  inclination  class  value  at  a  site  on  the  surface  will  be  the  site 
value  Xjl.  The  inclination  values  are  defined  by  the  set  X  of  three  integers 

X={  x  :  0 , 1  ,  2  }  (16) 

where  0  is  null,  1  principal  and  2  diagonal  local  inclinations. 

A  finite  spatial  neighborhood  of  radius  of  r=l  was  selected  to  investigate  the 
transition  of  local  inclination  class  xjl  to  xjt+l.  Therefore,  the  neighborhood 
configuration  number,  according  to  Equation  10, 

Vi1  =  22  x i_ i 1  +  21  x*  +  2°  xi+1[  <17) 


A  total  of  three  sites  compose  each  site  neighborhood  at  time  t.  From  Equation 
(11)  a  total  of  27  neighborhood  configurations  are  possible.  The  set  of 
probabilities  py<x  governing  the  transition  of  site  x;1  to  Xj,+  1  will  be  examined 
using  the  video  images  and  the  form  of  the  stochastic  rule  P,  which  is  now 

Xjt+1  =  P[  xj.i1  ,  Xj‘  ,  xj  +  i1  ;  py>o  ,  py,i  ,  py,2  J  (18) 
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*il+1  =  P(  yil  ;  Py,0  •  Py.l  •  Py,2  ] 


(19) 


o  r 


The  probabilities  for  each  neighborhood  y  satisfy  the  requirement  that 

Py.O  +  Py,l  +  Py,2  =  1  (20) 

As  discussed  above,  the  definition  of  three  inclination  values  does  not 
uniquely  specify  the  fracture  surface  topography.  Rather,  it  defines  and 
distinguishes  between  the  classes  of  inclinations  on  the  surface.  It  should  be 
noted  that  if  construction  of  the  actual  surface  topography  was  desired,  this 
could  be  accomplished  by  other  methods. 


Image  processing  algorithm 

The  set  Q  of  eight  images  was  obtained  for  processing  with  respect  to  the  set  C 
of  calibration  curves.  Each  of  the  eight  images  were  obtained  using  a 
different  light  source  from  the  set  L.  The  fracture  surface  plane  F  was  normal 
to  the  optical  axis  O  for  all  eight  fracture  surface  images.  Reduction  of  the 
eight  images  to  the  processed  image  of  the  surface  was  based  on  a  three  step 
procedure  to  identify  the  local  surface  inclination  class  at  each  pixel  from  the 
null,  principal  and  diagonal  inclinations  defined  previously  (Figure  7). 

The  first  step  of  the  procedure  used  the  calibration  curves  to  transform  the 
pixel  intensities  qtl  from  values  ranging  from  0  to  255  to  a  value  of  0  or  1  with 
respect  to  a  threshold  inclination  angle  vih  selected  by  the  user.  A  set  of 
threshold  angles  was  chosen  to  uniformly  span  a  range  from  zero  to  null 
saturation.  For  this  investigation,  the  range  was  from  0°  to  20°  at  4° 
increments.  The  threshold  inclination  lias  an  associated  threshold  intensity 
associated  with  each  calibration  curve  (Figure  8).  As  discussed  above,  six 
calibration  curves  for  each  of  the  two  specimen  rotation  axes  exhibited 
sensitivity  to  specimen  rotation.  These  twelve  curves  were  paired  with  the 
appropriate  fracture  surface  image  according  to  illumination  direction. 
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Schematic  of  pise!  intensity  to  local  surface  inclination  transformation 
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of  nonlinear  intensity  versus  specimen  orientation  data  to  binary  value 


Higher  pixel  intensity  indicated  an  orientation  toward  the  illumination  source. 
The  shape  of  the  calibration  curve  was  sufficiently  nonlinear  to  identify 
orientation  away  from  the  illumination  source  which  exceeded  the  threshold 
inclination.  Decreasing  sensitivity  of  the  calibration  curve,  corresponding  to 
increased  inclination  in  the  direction  of  illumination,  did  not  support 


discrimination  between  inclination  smaller  than  the  threshold  inclination 

and  inclination  greater  than  the  threshold  in  the  direction  of  the 

illumination. 

The  intensity  of  each  pixel  in  each  of  the  eight  images  was 
associated  calibration  curve  intensity,  consistent  with 
inclination  away  from  the  illumination  source,  i.e. 

compared  with  its 

its  sensitivity  to 

otil,e$  =  (X  [  qil,$  ,  cit,e0vth  1 

(20) 

9  =  0°  .90° 

(21) 

<t>  =  e  -  45° ,  e ,  e  +  45° . 

(0+180°  )-45°  ,  (9+180°  ) ,  (0+18O°)+45° 

(22) 

If  the  pixel  intensity  was  greater  than  the  threshold  intensity,  the  pixel  was 
given  a  transformed  intensity  value  of  0.  If  the  pixel  intensity  was  less  than 
the  threshold  intensity,  a  transformed  value  of  1  was  assigned.  A  value  of  0 

indicated  no  apparent  inclination  of  the  pixel  field  with  respect  to  the 

threshold.  A  value  of  1  indicated  significant  inclination  with  respect  to  the 

threshold  value.  Through  this  process,  each  portion  of  the  region  of  interest 
on  the  fracture  surface  associated  with  a  pixel  was  characterized  by  twelve 
binary  values  grouped  into  four  three-bit  words. 

The  second  step  in  the  image  processing  sequence  examined  the  ajl,ea  binary 
values  to  assess  whether  these  attributes  are  consistent  with  the  four  possible 
principal  inclinations.  These  binary  values  were  considered  in  sets  of  three 
adjacent  ajlte$  values,  as  summarized  in  Figure  9,  for  each  principal  direction 

which  would  be  considered.  A  concensus  for  a  candidate  principal  direction 
was  defined  according  to  a  majority  rule.  Each  of  the  four  three-bit  binary 
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words  were  transformed  to  a  single  four-bit  binary  value  in  this  manner 
(Table  1),  i.e. 


Pi'-e+di  =  (3  [  ait.e  0-450  ,  a^e  0 ,  ajt.e  0+450  ] 

(23) 

0  =  0°  ,  90° 

(24) 

<t»  =  e ,  &+18O0 

(25) 

Each  binary  value  was  associated  with  one  of  the  four  possible  principal 
directions. 

The  third  image  processing  step  examined  the  symmetries  of  the  four  possible 
principal  inclinations  in  order  to  arrive  at  consistent  representations  of  the 
null,  principal  and  diagonal  inclination  classes  from  these  values.  The 
assignment  of  an  inclination  class  to  each  pixel  was  based  on  the  rotational 
asymmetry  of  the  Pil,e+$  indicators  (Figure  10),  i.e. 

xjl  =  y[  b^oo  ,  bil,90o  ,  bj^igoo  ,  bilt270o  ]  (26) 

Four  null  values  suggest  a  null  inclination  class.  One  prinipal  value  and  three 
null  values  suggests  a  principal  inclination  class.  Two  principal  values  and 
two  null  values  suggest  a  diagonal  inclination  class  if  asymmetrically 
distributed.  Two  principal  values,  if  symmetrically  distributed,  offer  no 
concensus  and  therefore  suggest  a  null  inclination.  Three  out  of  four  or  four 
out  of  four  principal  inclination  values  are  inconsistent  with  a  well  defined 
principal  inclination  and  are  taken  to  be  null  inclination.  In  this  way,  the 
four-bit  binary  word  argument  is  transformed  to  a  1-bit  ternary  word 
indicating  inclination  class  (Table  2). 


Discussion  of  Processed  Images 

For  purposes  of  presentation,  the  null,  principal  and  diagonal  inclination 
values  were  represented  by  red,  green  and  blue  colors,  respectively.  Six 
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Table  1 


ail,04 

Pil,8+$ 

000 

0 

001 

0 

010 

0 

Oil 

1 

100 

0 

101 

1 

110 

1 

111 

1 

—  Transformation  of  3-bit  binary  inclination  words  to  1-bit  binary  candidate 
principal  inclination  value  by  majority  rule 
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x j *  =0  (NULL) 

(a) 

Fig.  10  —  Transformation  of  candidate  principal  inclinations  to  (a)  null  inclinations 
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X|<  =  1  (PRINCIPAL) 


o 

o 

o 

o 

Xj'  =  2  (DIAGONAL) 


(b)  (c) 

Fig.  10  —  (Continued)  Transformation  of  candidate  principal  inclinations  to  (b)  principal  and 

(c)  diagonal  inclinations 


30 


Xil 

0000 

0 

0001 

1 

0010 

1 

0011 

2 

0100 

1 

0101 

0 

0110 

2 

0111 

0 

1000 

1 

1001 

2 

1010 

0 

1011 

0 

1100 

2 

1101 

0 

1110 

0 

mi 

0 

Table  2  —  Transformation  of  4-bit  binary  candidate  principal  inclination  word  to 
ternary  inclination  value  from  symmetry  considerations 
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processed  images  were  generated  for  threshold  angles  yth  values  from  0°  to  20° 
at  4°  increments.  The  processed  images  are  displayed  in  the  bottom  row  of 
Figure  11. 

Threshold  angles  near  0°  presented  an  ambiguous  state  to  the  algorithm.  The 

threshold  decision  in  the  image  processing  algorithm  would  theoretically  be 
based  on  identical  or  nearly  identical  intensity  values  for  each  pixel  since  at 
0°  threshold  all  such  images  are,  in  principle,  degenerate.  Noise  in  the 
calibration  curves  dominated  the  threshold  decision  at  small  threshold  angles. 

As  the  threshold  angle  increased,  geometric  structure  emerged  on  the 
processed  images.  The  processed  images  demonstrated  a  complex  geometric 
grouping  of  local  inclination  classes  on  the  fracture  surface.  As  the  threshold 
inclination  angle  was  increased,  the  processed  image  tended  to  a  uniform  null 

inclination.  The  image  becomes  increasingly  saturated  by  the  red  color 
assigned  to  pixels  of  null  inclination.  This  was  consistent  with  the 
classification  of  each  pixel  site  as  a  null  inclination  with  respect  to  the 
reference  plane  R.  The  null  inclination,  therefore,  functions  as  a  reference 
and  background  upon  which  principal  and  diagonal  inclinations  are 
superimposed. 

The  dominant  trend  at  the  intermediate  values  of  threshold  inclination  angle 
(over  the  range  considered)  is  for  principal  inclinations  ,  indicated  by  green, 
to  be  grouped  together  as  "islands"  on  the  null  inclination  "sea"  of  red.  The 

blue  diagonal  inclination  class  never  appears  as  isolated  structures  on  the  null 
background.  All  blue  features  exist  as  aggregate  forms  on  the  green  principal 
inclination  structures  in  the  processed  image.  In  fact,  relatively  few  instances 
of  blue-red  boundaries  exist.  This  implies  that  the  diagonal  inclinations  are 
coupled  to  the  principal  inclinations  in  a  secondary  role  for  this  particular 
fracture  surface. 

Whether  the  diagonal  inclinations  are  a  cause,  an  effect  or  a  combined 
influence  on  the  formation  of  the  fracture  surface  cannot  be  inferred  from 
the  basis  of  this  data.  The  relative  proportion  and  relationship  of  diagonal 
inclination  areas  to  principal  inclination  areas  suggest  that  the  diagonal 
inclination  sites  are  dominated  by  smaller  scale  features  which  disrupt,  in 
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He.  1 1  —  The  processed  iracture  surlace  image,  a  pro; 
an.loni  <eed  synthetic  processed  image  and  a  random  nn 


some  maimer,  the  principal  inclination  features  on  the  fracture  surface  run 
parallel  to  the  direction  of  crack  growth.  The  diagonal  inclination  sites  would 
appear  to  be  spatially  infrequent  but  capable  of  strong  local  action.  Closer 

examination  of  these  sites  by  higher  resolution  fractographic  methods  could 
provide  additional  insight  into  the  relationship  between  microstructural  and 
mezzostructural  features. 

The  stochastic  cellular  automaton  rule  for  each  value  of  the  threshold  angle 
was  constructed  by  tabulating  the  xit+1  as  outputs  for  each  yjl  neighborhood 
configuration  as  input.  The  three  tabulated  outputs  for  each  neighborhood 
configuration  were  normalized.  Therefore,  the  probabilities  in  each  row 
summed  to  unity,  according  to  Equation  (IS).  Only  two  of  the  probabilities  in 
each  row  were  independent.  Each  resulting  rule  table  featured  27  rows,  one 
for  each  possible  input  configuration,  and  3  columns,  one  for  each  possible 
output.  Each  value  in  the  table  was  the  probability  that  the  input 
configuration  of  the  row  would  produce  the  output  configuration  of  the 
column.  Each  table  contained  54  independent  parameters. 

Of  the  27  neighborhood  configurations,  nine  are  symmetric  under  spatial 
reflection.  The  remaining  eighteen  can  be  divided  into  nine  pairs  of 

asymmetric  configurations  which  map  into  one  another  under  spatial 
reflection.  The  physical  nature  of  the  fracture  surface,  as  captured  by  the 
stochastic  cellular  automaton,  was  examined  by  performing  a  two  factor 

analysis  of  variance  on  the  three  output  occurrences  associated  with  each 
asymmetric  neighborhood  configurations  pair.  A  two  row  by  three  column 
table  was  created  for  each  of  the  nine  neighborhood  configuration  pairs.  The 

rows  corresponded  to  the  input  configurations  and  the  columns  to  the  output 
site  values.  Symmetry  conditions  suggest  that  the  output  probabilities 
associated  with  001  and  100  input  configurations,  for  example,  should  be  equal. 

F  factors  for  the  row  variance  to  random  variance  ratio,  obtained  for  each 

neighborhood  pair  at  each  value  of  the  threshold  angle,  were  in  42  of  45  cases 
insufficient  to  reject  a  null  hypothesis  of  symmetric  behavior  at  a  0.05 
significance  level  (Table  3).  In  one  case,  at  high  threshold  angle,  the  tabulated 

data  would  not  support  an  analysis  of  variance,  obtained  from  the  processed 

images  These  results  suggest  that  the  aggregate  local  physical  symmetry  of  the 
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Table  3  —  Two  factor  analysis  of  variance  F  values  for  row  variance  to  random  variance  ratio. 
Significance  level  of  0.05,  indicating  departure  from  translational  symmetry  with  respect  to 
neighborhood  configuration,  corresponds  to  F  value  of  18.5. 


35 


surface  is  supported  by  the  data  from  the  processed  images  and  the  stochastic 
cellular  automata  formalism. 


Generation  of  synthetic  images 

If  the  stochastic  cellular  automata  formalism  successfully  captures  the  small 
scale  interactions  and  large  scale  structure  of  the  fracture  surface,  iterative 
application  of  the  automaton  should  be  able  to  generate  synthetic  processed 
images.  Different  initial  conditions,  in  terms  of  site  values  at  zero  time  and 
random  number  sequences  implementing  the  rule  probabilities,  were 
considered  in  order  to  examine  the  robustness  of  the  stochastic  cellular 

automata.  Essential  features  of  the  processed  image,  such  as  the  relative 
proportions  of  null,  principal  and  diagonal  inclination  classes,  the  self 

organization  between  inclination  classes  and  the  connectivity  between 
inclination  class  regions  all  indicate  the  ability  of  the  stochastic  cellular 

automata  formalism  to  efficiently  and  compactly  capture  the  essence  of  the 

surface  topology. 

Two  types  of  synthetic  processed  images  were  generated.  In  the  first  type  of 
synthetic  processed  image,  the  zero  time  row  of  the  processed  physical  image 
was  used  to  seed  the  zero  time  row  of  the  synthetic  image.  Random  number 
sequences  were  used  to  support  the  rule  table  probability  decisions  required 
by  the  automaton.  For  the  second  type  of  processed  image,  the  zero  time  row  of 
the  synthetic  image  was  seeded  by  random  initial  values  in  the  same 
proportions  exhibited  by  the  processed  image.  Random  number  sequences 

were  again  used  to  support  the  required  rule  probability  decisions.  The 

synthetic  processed  images  are  shown  in  the  second  and  third  rows  of  Figure 

11.  The  self  organization  and  similarity  of  the  aggregate  structures  in  the 
processed  and  synthetic  images  are  apparent.  (One  synthetic  image  is  shown 
for  each  case  in  Figure  11.  Five  images  for  each  case,  further  illustrating  the 
similarity  of  image  structures,  are  shown  in  the  Appendix.) 

A  set  of  red.  green  and  blue  random  images  is  also  shown  in  the  top  row  of 

Figure  1 1  for  reference  and  comparison.  The  random  images  color  proportions 
were  selected  to  be  equal  to  the  color  proportions  of  the  respective  processed 
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image.  The  departure  of  the  processed  and  synthetic  processed  images  from 
randomness  is  evident. 


Concluding  remarks 

Experimental  measurement  of  local  fracture  surface  inclination  were  obtained 
from  a  set  of  video  images  for  intensity  versus  orientation  calibration  and  a  set 
of  video  images  for  processing.  The  images  for  processing  were  mapped  to  a  set 
of  inclination  classes  consistent  with  null,  principal  and  diagonal  inclinations 
defined  with  respect  to  the  macroscopic  surface  normal  and  crack  growth 
direction. 

The  global  character  of  the  surface,  in  terms  of  local  interactions  between 
null,  principal  and  diagonal  inclination  classes,  was  examined  through 
identification  of  stochastic  cellular  automaton  parameters  from  the  processed 
image  of  local  inclination  classes.  The  automaton  output  probabilities 
corresponding  to  input  neighborhood  configurations  equivalent  under  spatial 
reflection  were  found  to  be  similar.  Physical  consistency  of  the  parameters 
with  respect  to  symmetry  under  spatial  reflection  was  therefore  demonstrated. 

Synthetic  processed  images  were  generated  using  two  different  types  of  initial 
conditions.  The  first  set  of  synthetic  processed  images,  referred  to  as  processed 
seed  synthetic  processed  images,  were  generated  from  the  processed  local 
inclination  values  at  the  zero  time  sites.  The  image  boundary  corresponded  to 
the  direction  in  which  crack  propagation  generated  the  surface  in  the  image 
field.  The  stochastic  cellular  automaton  probabilities  were  invoked  with  a 
random  number  sequence  to  arrive  at  the  output  decision  for  each  input  local 
neighborhood  configuration.  A  total  of  five  different  processed  seed  synthetic 
images  were  generated. 

The  second  set  of  synthetic  processed  images,  referred  to  as  random  seed 
synthetic  processed  images,  were  generated  from  zero  time  sites  using  random 
seeding  in  proportion  to  the  processed  image  local  inclination  class 
probabilities.  Output  decisions  for  these  synthetic  images  were  also  obtained 
using  random  number  sequences.  Five  random  seed  synthetic  images  were 
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generated  from  five  different  sets  of  initial  conditions  and  five  different 
random  number  sequences. 

The  processed  and  synthetic  processed  images  exhibited  features  similar  to 
those  of  the  processed  image  obtained  from  experimental  data.  All  three  types 
of  images  show  considerable  self  organization  of  the  local  inclination  classes 

when  compared  to  random  images. 

The  processed  and  synthetic  processed  images  demonstrate  the  size  and  shape 

of  local  inclination  class  aggregation  on  the  fracture  surface.  The  fracture 
surface  is  a  manifestation  of  the  underlying  antecedant  damage  processes  at  a 
smaller  scale  in  the  material.  The  information  transport  in  the  stochastic 
cellular  automaton  reflects  a  balance  between  the  strength  of  the  small  scale 

local  interactions,  defined  by  the  input  neighborhood,  and  the  persistence  of 
large  scale  damage  structures  which  have  evolved  from  these  interactions 

Metallic  fracture  surfaces  are  complex  geometric  structures.  The  spatial  scales 
defining  material  damage  and  fracture  are  evolving  scales.  Both  are 

simultaneously  united  at  the  instant  of  surface  generation.  The  material 
damage  evolution  is  frozen  by  the  very  fracture  event  which  it  precipitates 
while  the  fracture  surface  is  defined  by  the  aggregate  state  of  the  damage  at 
fracture.  The  ability  to  identify  the  dominant  local  interaction  scale  and  the 
global  scale  of  self-organization  from  experiment  should  complement 
analytical  modeling  of  the  fracture  process.  This  method  of  analysis  serves  to 
identify  extraordinary  features  on  an  already  complex  background  for 
additional  scrutiny  by  higher  resolution  experimental  methods. 
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Appendix 


Tabulated  data  from  each  processed  image,  for  each  threshold  angle  value 

considered,  is  provided  in  Table  Al.  These  tables  were  used  to  determine  the 

parameters  for  the  stochastic  cellular  automata.  The  stochastic  cellular 
automata  rules  was  constructed  by  normalizing  the  x,t+1  as  outputs  with 

respect  to  the  yjl  neighborhood  configuration  as  input. 

As  discussed  in  the  main  text,  two  types  of  synthetic  processed  images  were 
generated.  In  the  first  type  of  synthetic  processed  image,  the  zero  time  row  of 
the  processed  physical  image  was  used  to  seed  the  zero  time  row  of  the 

synthetic  image.  For  the  second  type  of  processed  image,  the  zero  time  row  of 
the  synthetic  image  was  seeded  by  random  initial  values  in  the  same 

proportions  exhibited  by  the  processed  image.  Random  number  sequences 
were  used  to  suppon  the  rule  table  probability  decisions  in  both  cases.  Five 

synthetic  images  of  each  type,  based  on  the  stochastic  rule  parameters  derived 
from  the  processed  images  for  each  value  of  the  threshold  angle,  are  shown  in 
Figure  Al.  The  self  organization  and  similarity  of  the  aggregate  structures  in 

the  processed  and  synthetic  images  were  found  to  be  relatively  insensitive  to 

the  random  number  seed  or  to  the  initial  conditions  used. 
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Table  A1  —  (Continued)  Tabulated  output  site  values  with  respect  to  input  neighborhood 

configurations  for  threshold  angles  (ii)  4° 
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Table  A1  —  (Continued)  Tabulated  output  site  values  with  respect  to  input  neighborhood 
configurations  for  threshold  angles  of  (iii)  8° 
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Table  A1  —  (Continued)  Tabulated  output  site  values  with  respect  to  input  neighborhood 
configurations  for  threshold  angles  of  (v)  16° 
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